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Abstract 

We deal with the asymptotic enumeration of combinatorial structures on planar maps. 
Prominent instances of such problems are the enumeration of spanning trees, bipartite 
perfect matchings, and ice models. The notion of orientations with out-degrees prescribed 
by a function a : V — > N unifies many different combinatorial structures, including the 
afore mentioned. We call these orientations a-orientations. The main focus of this pa- 
per are bounds for the maximum number of a-orientations that a planar map with n 
vertices can have, for different instances of a. We give examples of triangulations with 
2.37" Schnyder woods, 3-connected planar maps with 3.209™ Schnyder woods and inner 
triangulations with 2.91™ bipolar orientations. These lower bounds are accompanied by 
upper bounds of 3.56™, 8™ and 3.97™ respectively. We also show that for any planar map 
M and any a the number of a-orientations is bounded from above by 3.73™ and describe 
a family of maps which have at least 2.598™ a-orientations. 

AMS Math Subject Classification: 05A16, 05C20, 05C30 

1 Introduction 

A planar map is a planar graph together with a crossing-free drawing in the plane. Many dif- 
ferent structures on planar maps have attracted the attention of researchers. Among them are 
spanning trees, bipartite perfect matchings (or more generally bipartite /-factors), Eulerian 
orientations, Schnyder woods, bipolar orientations and 2-orientations of quadrangulations. 
The concept of orientations with prescribed out-degrees is a quite general one. Remarkably, 
all the above structures can be encoded as orientations with prescribed out-degrees. Let a 
planar map M with vertex set V and a function a : V — > N be given. An orientation X of the 
edges of M is an a-orientation if every vertex v has out-degree a(v). For the sake of brevity, 
we refer to orientations with prescribed out-degrees simply as a-orientations in this paper. 

For some of the above mentioned structures it is not obvious how to encode them as ot- 
orientations. For Schnyder woods on triangulations the encoding by 3-orientations goes back 
to de Fraysseix and de Mendez [TO]. For bipolar orientations an encoding was proposed by 
Woods [40j and independently by Tamassia and Tollis [35j. Bipolar orientations of M are 
one of the structures which cannot be encoded as a-orientations on M, an auxiliary map 
M' (the angle graph of M) has to be used instead. For Schnyder woods on 3-connected 



*The conference version of this paper is to appear in the proceedings of WG'07 (LNCS) under the title "On 
the Number of a-Orientations" . 
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planar maps as well as bipartite /-factors and spanning trees Felsner p3] describes encodings 
as a-orientations. He also proves that the set of a-orientations of a planar map M can 
always be endowed with the structure of a distributive lattice. This structure on the set of 
a-orientations found applications in drawing algorithms in [3], [17] . and for enumeration and 
random sampling of graphs in [T3] . 

Given the existence of a combinatorial structure on a class M n of planar maps with n 
vertices, one of the questions of interest is how many such structures there are for a given 
map M £ M n - Especially, one is interested in the minimum and maximum that this number 
attains on the maps from M n . This question has been treated quite successfully for spanning 
trees and bipartite perfect matchings. For spanning trees the Kirchhoff Matrix Tree Theorem 
allows to bound the maximum number of spanning trees of a planar graph with n vertices 
between 5.02 n and 5.34 n , see [31] 128]. Pfaffian orientations can be used to efficiently calculate 
the number of bipartite perfect matchings in the planar case, see for example [23]. Kasteleyn 
has shown, that the k x £ square grid has asymptotically e°' 29 ' ki ~ 1.34 ke perfect matchings. 
The number of Eulerian orientations is studied in statistical physics under the name of ice 
models, see [2] for an overview. In particular Lieb [22J has shown that the square grid on the 
torus has asymptotically (8\/ r 3/9) ke ~ 1.53 Eulerian orientations and Baxter pQ has worked 
out the asymptotics for the triangular grid on the torus as (3y/3/2) M 2.598 M . 

In many cases it is relatively easy to see which maps in a class M n carry a unique object of 
a certain type, while the question about the maximum number is rather intricate. Therefore, 
we focus on finding the asymptotics for the maximum number of a-orientations that a map 
from M n can carry. The next table gives an overview of the results of this paper for different 
instances of M n and a. The entry c in the "Upper Bound" column is to be read as 0(c n ), in 
the "Lower Bound" column as Vl(c n ) and for the "~ c" entries the asymptotics are known. 



Graph class and orientation type 


Lower bound 


Upper bound 


a-orientations on planar maps 


2.598 


3.73 


Eulerian orientations 


2.598 


3.73 


Schnyder woods on triangulations 


2.37 


3.56 


Schnyder woods on the square grid 


» 3.209 


Schnyder woods on 3-connected planar maps 


3.209 


8 


2-orientations on quadrangulations 


1.53 


1.91 


bipolar orientations on stacked triangulations 


2 


bipolar orientations on outerplanar maps 


» 1.618 


bipolar orientations on the square grid 


2.18 


2.62 


bipolar orientations on planar maps 


2.91 


3.97 



The paper is organized as follows. In Section [2] we treat the most general case, where M n 
is the class of all planar maps with n vertices and a can be any integer valued function. We 
prove an upper bound which applies for every map and every a. In Section [2.31 we deal with 
Eulerian orientations. In Section [3.11 we consider Schnyder woods on plane triangulations and 
in Section 13.21 the more general case of Schnyder woods on 3-connected planar maps. We 
split the treatment of Schnyder woods because the more direct encoding of Schnyder woods 
on triangulations as a-orientations yields stronger bounds. In Section 13.21 we also discuss 
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the asymptotic number of Schnyder woods on the square grid. Section d] is dedicated to 
2-orientations of quadrangulations. In Section [5l we study bipolar orientations on the square 
grid, stacked triangulations, outerplanar maps and planar maps. The upper bound for planar 
maps relies on a new encoding of bipolar orientations of inner triangulations. In Section 16.11 
we discuss the complexity of counting a-orientations. In Section 16.21 we show how counting 
a-orientations can be reduced to counting (not necessarily planar) bipartite perfect matchings 
and the consequences of this connection are discussed as well. We conclude with some open 
problems. 



2 Counting a- Orient at ions 



A planar map M is a simple planar graph G together with a fixed crossing-free embedding planar map 
of G in the Euclidean plane. In particular, M has a designated outer (unbounded) face. We 
denote the sets of vertices, edges and faces of a given planar map by V, E, and J 7 , and their 
respective cardinalities by n, m and /. The degree of a vertex v will be denoted by d{v). 

Let M be a planar map and a : V — > N. An orientation X of the edges of M is an 
a-orientation if for all v G V exactly a{v) edges are directed away from v in X. a- 

Let X be an a-orientation of M and let C be a directed cycle in X. Define X c as the orientation 
orientation obtained from X by reversing all edges of C . Since the reversal of a directed 
cycle does not affect out-degrees the orientation X c is also an a-orientation of M. The plane 
embedding of M allows us to classify a directed simple cycle as clockwise (cw-cycle) if the cw-cycle 
interior, Int(C), is to the right of C or as counterclockwise (ccw-cycle) if Int(C) is to the left ccw-cycle 
of C. If C is a ccw-cycle of X then we say that X c is left of X and X is right of X c . Felsner left of 
proved the following theorem in [TJ]. right of 

Theorem 1 Let M be a planar map and a : V —* N. The set of a-orientations of M endowed 
with the transitive closure of the 'left of relation is a distributive lattice. 

The following observation is easy, but useful. Let M and a : V — > N be given, W C V and 
Ey/ the edges of M with one endpoint in W and the other endpoint in V \ W. Suppose all 
edges of Ey/ are directed away from W in some a-orientation Xq of M. The demand of W 
for YlweW a ( w ) outgoing edges forces all edges in Ey/ to be directed away from W in every 
a-orientation of M. Such an edge with the same direction in every a-orientation is a rigid 
edge. rigid edge 

We denote the number of a-orientations of M by r a (M). Let A4 be a family of pairs 
(M, a) of a planar map and an out-degree function. Most of this paper is concerned with 
lower and upper bounds for max( Mtt ) ejW r aM (M) for some family M.. In Section [2TT| we deal 
with bounds which apply to all M and a, while later sections will be concerned with special 
instances. 



2.1 An Upper Bound for the Number ^-Orientations 

A trivial upper bound for the number of a-orientations on M is 2 m as any edge can be directed 
in two ways. The following easy but useful lemma improves the trivial bound. 

Lemma 1 Let M be a planar map, A C E a cycle free subset of edges of M , and a a function 
a : V — > N. Then, there are at most 2 m_ l j4 l a-orientations of M. Furthermore, M has less 
than 4 n a-orientations. 
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Proof. Let X be an arbitrary but fixed orientation out of the 2 m ~' j4 ' orientations of the edges 
of E \ A. It suffices to show that X can be extended to an a-orientation of M in at most 
one way. We proceed by induction on \A\. The base case \ A\ = is trivial. If \A\ > 0, then, 
as A is cycle free, there is a vertex v, which is incident to exactly one edge e G A. If v has 
out-degree a{v) respectively a{v ) — 1 in X, then e must be directed towards respectively away 
from v. In either case the direction of e is determined by X, and by induction there is at 
most one way to extend the resulting orientation of E \ [A — e) to an a-orientation of M. 
If v does not have out-degree a(v ) or a(v) — 1 in X, then there is no extension of X to an 
a-orientation of M. The bound 2 m ~ n+1 < 4 n follows by choosing A to be a spanning forest 
and applying Euler's formula. □ 

A better upper bound for general M and a will be given in Proposition [TJ The following 
lemma is needed for the proof. 

Lemma 2 Let M be a planar map with n vertices that has an independent set I2 ofn<i vertices 
which have degree 2 in M . Then, M has at most (3n — 6) — [ri2 — 1) edges. 

Proof. Consider a triangulation T extending M and let B be the set of additional edges, i.e., 
of edges of T which are not in M. If n = 3 the conclusion of the lemma is true and we may thus 
assume n > 3 for the rest of the proof. Hence, there are no vertices of degree 2 in T, and every 
vertex of I2 must be incident to at least one edge from B. If there is a vertex v G I2, which is 
incident to exactly one edge from B, then v and its incident edges can be deleted from I2, from 
M and from T, whereby the result follows by induction. The last case is that all vertices of I2 
have at least two incident edges in B. Since every edge in B is incident to at most two vertices 
from I 2 it follows that \I 2 \ < \B\. Therefore, \E M \ = \E T \ - \B\ < \E T \ - \I 2 \ = (3n-6) -n 2 . 

□ 

Remark. It can be seen from the above proof, that i^2,n 2 pl us the edge between the two 
vertices of degree 712 is the unique graph to which only n2 — 1 edges can be added. For every 
other graph at least 712 edges can be added. 

Proposition 1 Let M be a planar map, a : V — > N, and L = L\ U L2 an independent set of 
M, where L2 is the subset of degree 2 vertices in I . Then, M has at most 

^'■n(^0 « 

a- orientations. 

Proof. We may assume that M is connected. Let Mj, for i = 1, . . . c, be the components of 
M — L. We claim that M has at most (3n — 6) — (c — 1) — (I/2I — 1) edges. Note, that every 
component C of M — L must be connected to some other component C via a vertex v € / 
such that the edges vw and vw 1 with w G C and w' G C form an angle at v. As to and w' 
are in different connected components the edge ww' is not in M and we can add it without 
destroying planarity. We can add at least c — 1 edges not incident to / in this fashion. Thus, 
by Lemma [2] we have that m + (c— 1) < 3n — 6 — (L2 — 1). 

Let S' be a spanning forest of M — /, and let S be obtained from S' by adding one edge 
incident to every v £ L. Then, S is a forest with n — c edges. By Lemma Q] M has at most 
2m-\s\ Q,-orientations and by Lemma [2] 

m-\S\< (3n - 6) - (c - 1) - (|I 2 | - 1) - (n - c) = 2n - 4 - \L 2 \. 
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For every vertex v G I\ there are 2 d ^~ 1 possible orientations of the edges of M — S at v. 
Only the orientations with a{v ) or a(v) — 1 outgoing edges at v can potentially be completed 
to an ce-orientation of M. Since I\ is an independent set it follows that M has at most 

2-15, . -TT _J_ ((d(v) - 1\ (d(v) - 1\\ 2n _ 4 „ |/2| „ /^_(d(t,)\\ 
a-orientations. □ 



Corollary 1 Let M be a planar map and a : V — > N. Then, M has at most 3.73™ a- 
orientations. 

Proof. Since M is planar the Four Color Theorem implies, that it has an independent set / 
of size |/| > n/4. Let 1%, I2 be as above. Note, that for d(v) > 3 

d(v)\ 1 ( d(v) \ 3 



2 d(v)-i \a(v)J ~ 2 d W-! VL rf W/ 2 J/ ~ 4 
Thus, the result follows from Proposition dj as 

2 2„_4-|/ 2 | ^ ^Y^' < 2 2„-4 ( 3y ^ . } T;r 



□ 



Remark. The best lower bound for general a and M, which we can prove, comes from 
Eulerian orientations of the triangular grid, see Section [2.31 



2.2 Grid Graphs 

Enumeration and counting of different combinatorial structures on grid graphs have received 
a lot of attention in the literature, see e.g. [21 [71 [22]. In Section [2.31 we present a family of 
graphs that have asymptotically at least 2.598™ Eulerian orientations. This family is closely 
related to the grid graph, and throughout the paper we will use different relatives of the grid 
graph to obtain lower bounds. We collect the definitions of these related families here. 
The grid graph Gj- n with k rows and £ columns is defined as follows. The vertex set is 

V k ,i = {(i,j) \ l<i<k, l<j<£}. 

The edge set = Ej? e U E^ £ consists of horizontal edges 

= {{(*> I l<i<k, 1 < j < ^ - l} 

and vertical edges 

Elt = {{(i,j),(i + hj)} \ l<i<k-l, l<j <£}. 

We denote the ith vertex row by = | 1 < j < £} and the jth vertex column by 

Vf = I 1 < i < k}. The jth edge column Ej is defined as E^ = + 1)} | 

1 < i < k}. The number of bipolar orientations of Gke is studied in Section [5.11 
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The grid on the torus e is obtained from Gk+i/+i by identifying and (k + l,i) 
as well as (j, 1) and (j, £ + 1) for all i and j, see Figures Q] (a) and (b). Edges of the form 
{(i, 1), (i,£)} are called horizontal wrap-around edges while those of the form {(1, j), (k,j)} 
are the vertical wrap-around edges. Note that Gk,e can be obtained from G\ t by deleting the 
k horizontal and the £ vertical wrap-around edges. 

Lieb [22] shows that GTe has asymptotically (8\/3/9) fc ^ Eulerian orientations. His analysis 
involves the calculation of the dominant eigenvalue of a so-called transfer matrix, see also 
Section [H 




Figure 1: Two illustrations of Gj A , the augmented grid G44, and the quadrangulation G^ 4 . 

We consider the number of Schnyder woods on the augmented grid G* k « in Section 13.21 
see Figure Q] (c) . The augmented grid is obtained from Gf- £ by adding a triangle with vertices 
{01,02,03} to the outer face. The triangle is connected to the boundary vertices of the grid 
as follows. The vertex a\ is adjacent to all vertices of V^, 0,2 is adjacent the vertices from 
Vf and 03 to the vertices from V k R U V{ . 

When we consider 2-orientations in Section H] we use the quadrangulation G ki , see Fig- 
ured (d). It is obtained from the grid Gk,i by adding one vertex Vco to the outer face which 
is adjacent to every other vertex of the boundary such that (1, 1) is not adjacent to Voo. For 
k and £ even this graph is closely related to the torus grid G^ e , which can be obtained from 
G^ e by reassigning end vertices of edge as follows. 

{(i,jVoo}^{(i,j),(M} ?<j<£ {{k,j), Voo } ^{(k,j),(i,j)} 2<j<e 
{(i,i), Voo } ^{(i,i),(i,e)} 2<i<k {(i,e), Voo } ^{(i,e),(i,i)} 2<i<k 

Since k, £ are even this does not create parallel edges and the resulting graph is G\ t minus 
the edges ei ={(1,1), (1,^)} and e 2 = {(1, 1), (k, 1)}. 

We also use the triangular grid T^i in Sections 12.31 and 15.21 It is obtained from G^ £ 
by adding the diagonal edges (i — 1, j + 1)} for 2 < i < k and 1 < j < £ — 1, see 

Figure [2] (a) . The augmented triangular grid TZ t , which we need in Section [3TT1 is obtained 
in the same way from G* k e , see Figure [71 

The terms vertex row, vertex column and edge column are used for the triangular grid 
analogously to the definition above for Gp. g. 

We also use the triangular grid on the torus T k e , see Figure [2] (b) . We adopt the definition 
from [I], therefore it differs slightly from that of the square grid on the torus. More precisely, 
instead of identifying vertices (i,£ + 1) and (i, 1) we identify vertices (i, £ + 1) and (i — 1, 1) 
(and (1,^ + 1) with (k, 1)) to obtain T kl from Tf./- This boundary condition is called helical. 
The wrap-around edges are defined analogously to the square grid case. 

Baxter [T] was able to determine the exponential growth factor of Eulerian orientations 
of T ki as k,£ — > 00. Baxter's analysis uses similar techniques as Lieb's [22] and yields an 

asymptotic growth rate of (3\/3/2) fc ^. 



6 




(b) 



(3,1) (3,2) (3,3) (3,4) 




(3,1) (3,2) (3,3) (3,4) 



J 3,4 



Figure 2: The triangular grid T^, and T% 4 . 



2.3 A Lower Bound Using Eulerian Orientations 

Let M be a planar map such that every v € V has even degree and let a be defined as 
a(v) = d(v)/2, \/v £ V. The corresponding a-orientations of M are known as Eulerian 
orientations. Eulerian orientations are exactly the orientations which maximize the binomial 
coefficients in equation ([!]). The lower bound in the next theorem is the best lower bound we 



have for max 
made for a. 



(M,a)eM r «M 



(M), where A4 is the set of all planar maps and no restrictions are 



Eulerian 
orientations 



Theorem 2 Let M n denote the set of all planar maps with n vertices and £{M) the set of 
Eulerian orientations of M G M. n - Then, for n big enough, 



2.59™ < (3V3/2) 



kf 



< max \£{M)\ < 3.73". 

M&Mn 



Proof. The upper bound is the one from Corollary [TJ For the lower bound consider the trian- 
gular torus grid TT,. As mentioned above Baxter pQ was able to determine the exponential 
growth factor of Eulerian orientations of T?o as k,£ — > oo. Baxter's analysis uses eigenvector 
calculations and yields an asymptotic growth rate of (3y/3/2) ke . This graph can be made into 
a planar map Tut by introducing a new vertex which is incident to all the wrap-around 
edges. This way all crossings between wrap-around edges can be substituted by Voo. As 
every Eulerian orientation of yields a Eulerian orientation of T^~ £ this graph has at least 
(3a/3/2) h > 2.598 H Eulerian orientations for k,£ big enough. □ 



3 Counting Schnyder Woods 

Sclmyder woods for triangulations have been introduced as a tool for graph drawing and graph 
dimension theory in |32l 133] . Schnyder woods for 3-connected planar maps are introduced 
in [12]. Here we review the definition of Schnyder woods and explain how they are encoded 
as a-orientations. For a comprehensive introduction see e.g. |13| . 

Let M be a planar map with three vertices a\, 02, 03 occurring in clockwise order on the 
outer face of M. A suspension M a of M is obtained by attaching a half-edge that reaches 
into the outer face to each of these special vertices. special 

Let M a be a suspended 3-connected planar map. A Schnyder wood rooted at 01,02,03 is vertices 

Schnyder 

7 wood 



an orientation and coloring of the edges of M a with the colors 1,2,3 satisfying the following 
rules. 

(Wl) Every edge e is oriented in one direction or in two opposite directions. The directions 
of edges are colored such that if e is bidirected the two directions have distinct colors. 

(W2) The half-edge at a% is directed outwards and has color i. 

(W3) Every vertex v has out-degree one in each color. The edges ei, e2, &z leaving v in colors 
1,2,3 occur in clockwise order. Each edge entering v in color i enters v in the clockwise 
sector from ej + i to see Figure [3] (a). 

(W4) There is no interior face the boundary of which is a monochromatic directed cycle. 



Figure 3: The left part shows edge orientations and edge colors at a vertex, the right part 
two different Schnyder woods with the same underlying orientation. 

In the context of this paper the choice of the suspension vertices is not important and we 
refer to the Schnyder wood of a planar map, without specifying the suspension explicitly. 

Let M a be a planar map with a Schnyder wood. Let denote the digraph induced by 
the directed edges of color i. Every inner vertex has out-degree one in T\ and in fact Tj is a 
directed spanning tree of M with root 

In a Schnyder wood on a triangulation only the three outer edges are bidirected. This 
is because the three spanning trees have to cover all 3n — 6 edges of the triangulation and 
the edges of the outer triangle must be bidirected because of the rule of vertices. Theorem [3] 
says, that the edge orientations together with the colors of the special vertices are sufficient 
to encode a Schnyder wood on a triangulation, the edge colors can be deduced, for a proof 
see PS|. 

Theorem 3 Let T be a plane triangulation, with vertices a%, a,2, occuring in clockwise order 
on the outer face. Let oct{v) := 3 if v is an internal vertex and axifli) '■= for i = 1,2,3. 
Then, there is a bijection between the Schnyder woods of T and the ay- orientations of the 
inner edges ofT. 

In the sequel we refer to an a^-orientation simply as a 3-orientation. Schnyder woods on 
3-connected planar maps are in general not uniquely determined by the edge orientations, see 
Figure [3] (b). Nevertheless, there is a bijection between the Schnyder woods of a 3-connected 
planar map M and certain a-orientations on a related planar map M, see |14| . 

In order to describe the bijection precisely, we first define the suspension dualM* ofM a , suspension 
which is obtained from the dual M* of M as follows. Replace the vertex v^, which represents dual 
the unbounded face of M in M* , by a triangle on three new vertices bi, 62, 63. Let Pi be the 
path from aj-i to a^+i on the outer face of M which avoids dj. In M a the edges dual to 



(a) 




2 
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those on Pi are incident to bi instead of v^. Adding a ray to each of the bi yields M a . An 
example is given in Figured! 




Figure 4: A Schnyder wood, the primal and the dual graph, the oriented primal dual com- 
pletion and the dual Schnyder wood. 

Proposition 2 Let M a be a suspended planar map. There is a bijection between the Schnyder 
woods of M a and the Schnyder woods of the suspension dual M a . Figure [21 illustrates how 
the coloring and orientation of a pair of a primal and a dual edge are related. 




Figure 5: The three possible oriented colorings of a pair of a primal and a dual edge. 

The completion M of M a and M a * is obtained by superimposing the two graphs such 
that exactly the primal dual pairs of edges cross, see Figure HI In the completion M the 
common subdivision of each crossing pair of edges is replaced by a new edge-vertex. Note 
that the rays emanating from the three special vertices of M a cross the three edges of the 
triangle induced by &i, &2, 63 and thus produce edge vertices. The six rays emanating into the 
unbounded face of the completion end at a new vertex Voo placed in this unbounded face. A 
pair of corresponding Schnyder woods on M a and M a induces an orientation of M which is 
an as-orientation where 

{3 for primal and dual vertices 
1 for edge vertices 
for Voo. 

Note, that a pair of a primal and a dual edge always consists of a unidirected and a bidirected 
edge, which explains why as{v e ) = 1 is the right choice. Theorem [4] says, that the edge 
orientations of M are sufficient to encode a Schnyder wood of M a , the edge colors can be 
deduced, for a proof see [H]. 

Theorem 4 The Schnyder woods of a suspended planar map M a are in bijection with the 
as -orientations of M. 

In the rest of this section we give asymptotic bounds for the maximum number of Schnyder 
woods on planar triangulations and 3-connected planar maps. We treat these two classes 
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separately because the more direct bijection from Theorem [3] allows us to obtain a better 
upper bound for Schnyder woods on triangulations than for the general case. We also have a 
better lower bound for the general case of Schnyder woods on 3-connected planar maps than 
for the restriction to triangulations. 

Stacked triangulations are plane triangulations which can be obtained from a triangle Stacked tri- 
by iteratively adding vertices of degree 3 into bounded faces. The stacked triangulations angulations 
are exactly the plane triangulations which have a unique Schnyder wood and we have a 
generalization of this well-known result for general 3-connected planar maps, which we state 
here without a proof. 

Theorem 5 All 3-connected planar maps, which have a unique Schngder wood, can be con- 
structed from the unique Schngder wood on the triangle bg the six operations show in Figure 
read from left to right. 




Figure 6: Using the three primal operations in the first row and their duals in the second row 
every graph with a unique Schnyder wood can be constructed. 



3.1 Schnyder Woods on Triangulations 

Bonichon [3j found a bijection between Schnyder woods on triangulations with n vertices and 
pairs of non-crossing Dyck-paths, which implies that there are C n+ 2C n — C^ +1 Schnyder woods 
on triangulations with n vertices. By C n we denote the nth Catalan number C n = ( / 
Hence, asymptotically there are about 16 n Schnyder woods on triangulations with n vertices. 
Tutte's classic result [38J yields that there are asymptotically about 9.48 n plane triangulations 
on n vertices. See [27] for a proof of Tutte's formula using Schnyder woods. The two results 
together imply that a triangulation with n vertices has on average about 1.68™ Schnyder 
woods. The next theorem is concerned with the maximum number of Schnyder woods on a 
fixed triangulation. 

Theorem 6 Let T n denote the set of all plane triangulations with n vertices and S{T) the 
set of Schngder woods of T £ T n . Then, 

2.37™ < max|«S(T)| < 3.56™. 
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The upper bound follows from Proposition Q] by using that for d(v) > 3 

( d ( v )\ 2 i-d(v) < 5 
V 3 J " 8' 

For the proof of the lower bound we use the augmented triangular grid T£g. Figure [7J shows 
a canonical Schnyder wood on T£ t in which the vertical edges are directed up, the horizontal 
edges to the right and diagonal ones left-down. 




Figure 7: The graphs T| 5 with a canonical Schnyder wood and with the additional edges 
simulating Baxter's boundary conditions. 

Instead of working with the 3-orientations of T kl we use a* -orientations of T^/ where 

{3 if2<i<A;-land2<i<^-l 
l if (i,j)e{(l,l),(M),(M)} 
2 otherwise. 

For the sake of simplicity, we refer to a* -orientations of T^ e as 3-orientations. 

Intuitively, Tk i promises to be a good candidate for a lower bound because the canonical 
orientation shown in Figure [7J on the left has many directed cycles. We formalize this intuition 
in the next proposition, which we restrict to the case k = I only to make the notation easier. 

Proposition 3 The graph T£ , has at least 2 5 / 4 ( fc-1 ) 2 Schnyder woods. For k big enough T? k 
has 

2.37 k2+3 < |<S(T fe * fe )| < 2.599 fc2+3 

Proof. The face boundaries of the triangles of Tf. & can be partitioned into two classes C and 
C of directed cycles, such that each class has cardinality (k — l) 2 and no two cycles from the 
same class share an edge. Thus, a cycle C £ C shares an edge with three cycles from C if it 
does not share an edge with the outer face of T^k and otherwise it shares an edge with one 
or two cycles from C. 

For any subset D of C reversing all the cycles in D yields a 3-orientation of Tf. and we 
can encode this orientation as a 0-1-sequence of length (fc — l) 2 . After performing the flips of 
a given 0-1-sequence a, an inner cycle C € C is directed if and only if either all or none of 
the three cycles sharing an edge with C' have been reversed. If C G C is a boundary cycle, 
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then it is directed if and only none of the adjacent cycles from C has been reversed. Thus the 
number of different cycle flip sequences on C U C is bounded from below by 



ae{o,i} (fc - 1)2 

Here Xc (a) is an indicator function which takes value 1 if C is directed after performing the 
flips of a and otherwise. 

We now assume that every a £ {0, l}( fc_1 ) is chosen uniformly at random. The expected 
value of the above function is then 

Ef2^-^cl = ^ \^ 2^-' c " sC ' x c( a ) 

ae{o,i} (,i - 1)2 

Jensen's inequality E[</?(X)] > (/?(E[X]) holds for a random variable X and a convex 
function ip. Using this we derive that 

E[2^^c'] > 2 E E*cl = 2 E p [ c " fii PP ablc ]. 

The probability that C is flippable is at least 1/4. For C which does not include a boundary 
edge the probability depends only on the three cycles from C that share an edge with C and 
two out of the eight flip vectors for these three cycles make C flippable. A similar reasoning 
applies for C including a boundary edge. Altogether this yields that 

ae{o,i}( k -v 2 

Different cycle flip sequences yield different Schnyder woods. The orientation of an edge is 
easily determined. The edge direction is reversed with respect to the canonical orientation if 
and only if exactly one of the two cycles on which it lies has been flipped. We can tell a flip 
sequence apart from its complement by looking at the boundary edges. 

For the upper bound we use Baxter's result for Eulerian orientations on the torus 7jT« 
(see Sections 12.21 and I2.3[) . Every 3-orientation of T^i plus the wrap-around edges, oriented 
as shown in Figure [7] on the right, yields a Eulerian orientation of T^ t We deduce that T£g 
has at most 2.599 n Schnyder woods. □ 

Remark. Let us briefly come back to the number of Eulerian orientations of Tut, which 
was mentioned in Sections 12.21 and 12.31 and in the above proof. There are only 2 2 ( fc+ ^ _1 
different orientations of the wrap-around edges. By the pigeon hole principle there is an 
orientation of these edges which can be extended to a Eulerian orientation of in asymp- 
totically (3>/3/2)** ways. Thus, there are out-degree functions ctke for Tk e such that there are 
asymptotically 2.598* a^-orientations. Note, however, that directing all the wrap-around 
edges away from the vertex to which they are attached in Figure [7] induces a unique Eulerian 
orientation of T^g. 

We have not been able to specify orientations of the wrap-around edges, which allow to 
conclude that T^ g has (3V3/2) ki 3-orientations with these boundary conditions. In particular 
we have no proof that Baxter's result also gives a lower bound for the number of 3-orientations. 
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3.2 Schnyder Woods on the Grid and 3-Connected Planar Maps 

In this section we discuss bounds on the number of Schnyder woods on 3-connected pla- 
nar maps. The lower bound comes from the grid. The upper bound for this case is much 
larger than the one for triangulations. This is due to the encoding of Schnyder woods by 
3-orientations on the primal dual completion graph, which has more vertices. We summarize 
the results of this section in the following theorem. 

Theorem 7 Let be the set of 3-connected planar maps with n vertices and S{M) denote 
the set of Schnyder woods of M G M.^. Then, 

3.209" < max \S(M)\ < 8 n . 

MeMl 

The example used for the lower bound is the square grid Gk,e- 

Theorem 8 For k, I big enough the number of Schnyder woods of the augmented grid Gt g is 
asymptotically \S(G* k/ )\ m 3.209 H . 

Proof. The graph induced by the non-rigid edges in the primal dual completion G* k « of G* k » 
is G < 2fc_i i 2£-i — {2k — 1, 1). This is a square grid of roughly twice the size as the original 
and with the lower left corner removed. The rigid edges can be identified using the fact that 
c*s{voc) = and deleting them induces a' s on G^fe-i^-i — (2fe — 1,1). The new a' s only 
differs from as for vertices, which are incident to an outgoing rigid edge, and it turns out, 
that ct'c(y) = d{v) — 1 for all primal or dual vertices and a' s {v) = 1 for all edge vertices of 
G2k-l,2i—l ~ (2/s — 1, 1). Thus, a bijection between a^-orientations and perfect matchings 
of G2k-i,2k-£ ~ (2& — 1,1) is established by identifying matching edges with edges directed 
away from edge vertices. The closed form expression for the number of perfect matchings of 
G2k-i,2k-e — (2fc — 1, 1) is known (see [21]) to be 

nn(4-2cos^-2cos^). 

i=l j=i \ / 

The number of perfect matchings of Cr2fc— l 2^-1 — {2k — 1, 1) is sandwiched between that 
of G2k-2,2£-2 and that of G2k,2i- Therefore the asymptotic behavior is the same and in [23] . 
the limit of the number of perfect matchings of G^fci^, denoted as &(2k,2£), is calculated to 
be 

log$(2fc,2l) log2 1 rr 2 a/^j^noa 

fc ]fco 2k- 2£ = — + 4^ J I l0g(C0S {X) + C ° S {v))dxdV W °- 29 ' 

This implies that G k « has asymptotically e 4 ' 0,29 '^ ~ 3.209 fc ^ Schnyder woods. □ 

Remark. In [36] Temperley discovered a bijection between spanning trees of Gk,£ and perfect 
matchings of G^fc-i^-i — {2k — 1, 1). Thus, Schnyder woods of G* kl are in bijection with 
spanning trees of Gk£, see Figure El This bijection can be read off directly from the Schnyder 
wood: the unidirected edges not incident to a special vertex form exactly the related spanning 
tree. Encoding both, the Schnyder woods and the spanning trees, as a-orientations also gives 
an immediate proof of this bijection. 

We now turn to the proof of the upper bound stated in Theorem [71 The proof uses 
the upper bound for Schnyder woods on plane triangulations, see Theorem [6l We define a 
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Figure 8: A Schnyder wood on G44, the reduced primal dual completion G?j — (7,1) with 
the corresponding orientation and the associated spanning tree. 

triangulation Tm such that there is in injective mapping of the Schnyder woods of M to the 
Schnyder woods of Tm- The triangulation Tm is obtained from M by adding a vertex vp to 
every face F of M with \F\ > 4, see Figure [H The generic structure of a bounded face of a 
Schnyder wood is shown on the left in the top row of Figure [9l for a proof see |13j . The three 
edges, which do not lie on the boundary of the triangle, are the special edges of F. special edges 

3232 2323 23 




Figure 9: A Schnyder wood on a map M induces a Schnyder wood on Tm- The three special 
edges of a face are those, which do not lie on the black triangle. 

A vertex vp is adjacent to all the vertices of F. A Schnyder wood of M can be mapped 
to a Schnyder wood of Tm using the generic structure of the bounded faces as shown in 
Figure [9j The green-blue non-special edges of F become green unidirected. Their blue parts 
are substituted by unidirected blue edges pointing from their original start-vertex towards 
vp. Similarly the blue-red non-special edges become blue unidirected and the red-green ones 
red unidirected. Three of the edges incident to vp are still undirected at this point. They are 
directed away from vp and colored in accordance with (W3). 

Let two different Schnyder woods be given that have different directions or colors on an 
edge e. That the map is injective can be verified by comparing the edges on the boundary of 
the two triangles on which the edge e lies in Tm- 
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Thus, it suffices to bound the number of Schnyder woods of Tm- We do this by specializing 
Proposition [TJ We denote the set of vertices of Tm that correspond to faces of size 4 in M 
by F 4 and its size by and similarly F> 5 and /> 5 are defined. Note that / = F4 U F> 5 is an 
independent set and Tm has a spanning tree in which all the vertices from / are leaves. Let 
tit denote the number of vertices of Tm • Then, Tm has at most 

2— — n (^(f 0) 5 4 " +/4+/!i • (2-)" ■ (§)*" - 4 " • 2 " • (§)* (4) 

Schnyder woods. Note that n + f& + />5 + f& + 2/>5 < m + f& + 2/>5 < 3n — 6 which implies 
that /4 + |/>5 < n. Maximizing equation (4) under this condition yields that the maximum 
8 n is attained when f± = n. Thus M has no more than 8 n Schnyder woods. □ 

The proof of the lower bound 3.209 n involves the result about the number of perfect match- 
ings of the square grid. This result makes use of non-combinatorial methods. Therefore, we 
complement this bound with a result for another graph family, which uses a straight-forward 
analysis, but still yields that these graphs have more Schnyder woods than the triangular 
grid, see Section [3TTI 

The graph we consider is the filled hexagonal grid Hk,£, see Figure [101 Neglecting bound- filled 
ary effects he hexagonal grid has twice as many vertices as hexagons. This can be seen by hexagonal 
associating with every hexagon the vertices of its northwestern edge. Thus, neglecting bound- grid 
ary effects, the filled hexagonal grid has five vertices per hexagon. The boundary effects will 
not hurt our analysis because H^i has only 2{k + i) boundary vertices but 5 • k£ + 2(k + I) 
vertices in total. 

Proposition 4 For k, t big enough the filled hexagonal grid H/-^ has 

2.63 n < \S(H k>t )\ < 6.07 n . 




Figure 10: The filled hexagonal grid -£^3, a Schnyder wood on this grid and the primal-dual 
suspension of a hexagonal building block of H^g. Primal vertices are black, face vertices blue 
and edge vertices green. 

Proof. We count how many different orientations we can have on a filled hexagon. We do 
this using the bijection from Theorem [H The right part of Figure [10] shows a feasible as- 
orientation of a filled hexagon. Note that this orientation is feasible on the boundary when 
we glue together the filled hexagons to a grid and add a triangle of three special vertices 
around the grid. We flip only boundary edges of a hexagon which belong to a 4-face in this 
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hexagon. As these edges belong to a triangle in the hexagon on their other side, the cycle 
flips in any two filled hexagons can be performed independently. 

Let us now count how many orientations a filled hexagon admits, see the right part of 
Figure [10] for the definition of the cycles Ci, 62,63 and C4. If the 6-cycle induced by the 
central triangle is directed as shown in the rightmost part of Figure [TU] then we can flip either 
C\ or C2 and if C2 is flipped, C3 can be flipped as well. This yields 4 3 orientations, as the 
situation is the same at the other two 4-faces of the hexagon. If the 6-cycle is flipped the 
same calculation can be done with C3 replaced by C4. This makes a total of 2 • 4 3 = 128 
orientations per filled hexagon. That is, there are at least 128 > 2.639 5 ' fc '^ orientations of 

We start the proof of the upper bound by collecting some statistics about H^^. As 
mentioned above, Hkl has n = 5 • k ■ t interior vertices, 12 • kl edges and 7 • kt faces. Thus, 
the primal-dual completion has 48 • kl edges. There is no choice for the orientation of the 
edges incident to the 3 • 4/7 • / = 12 • ki face vertices of triangles. We can choose a spanning 
tree T on the remaining 5 • kl + 12 • kt + 3 • kt vertices such that all face vertices are leafs 
and proceed as in the proof of Proposition [IJ but using that we know the number of edges 
exactly. Since in the independent set of the remaining face vertices all of them have degree 
4 and required out-degree 3, they contribute a factor of 1/2 each. Thus, there are at most 
2(48-12-20)*^ . 2-3-W = 2 13-M < 6.07 n Schnyder woods on H k/ . □ 



4 Counting 2-Orientations 



Felsner et al. [16] present a theory of 2-orientations of plane quadr angulations, which shows 
many similarities with the theories of Schnyder woods for triangulations. A quadrangulation 
is a planar map such that all faces have cardinality four. A 2-orientation of a quadrangulation 
Q is an orientation of the edges such that all vertices but two non-adjacent ones on the outer 
face have out-degree 2. 

In [15] it is shown that 2-orientations on quadrangulations with n inner quadrangles are 
counted by the Baxter-number B n+ \. Hence asymptotically there are about 8 n 2-orientations 
on quadrangulations with n vertices. Tutte gave an explicit formula for rooted quadrangu- 
lations. A bijective proof of Tutte's formula is contained in the thesis of Fusy [18J. The 
formula implies that asymptotically there are about 6.75 n quadrangulations on n vertices. 
The two results together yield that a quadrangulation with n vertices has on average about 
1.19™ 2-orientations. 

We now give a lower bound for the number of 2-orientations of GrP*. The proof method 
via transfer matrices and eigenvalue estimates comes from Calkin and Wilf [7]. There it is 
used for asymptotic enumeration of independent sets of the grid graph. Let Z(Q) denote the 
set of all 2-orientations of a quadrangulation Q, with fixed sinks. 



2- 

orientation 
of a quad- 
rangulation 



Proposition 5 For k, t big enough has 



1.537" < \Z(GV e )\ < 



kl) 



>/3/9) M < 1.5397 M . 



Proof. We consider 2-orientations of GrP» with sinks (1, 1) and v^. These 2-orientations 
induce Eulerian orientations of GTp- The wrap-around edges inherit the direction of the 
respective edges incident to Voo (see Section [272]) and e±, ei are directed away from (1,1). 
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Figure 11: A 2-orientation of Ggg, the corresponding Eulerian orientation X of Gg'g and an 
alternating orientation of Gj 4 that can be extended to X. 



Therefore GtV has at most as many 2-orientations as G k « has Eulerian orientations, which 
implies the claimed upper bound. 

Conversely a Eulerian orientation of G ki in which the wrap-around edges have these 
prescribed orientations induces a 2-orientation of Gj?g. Such Eulerian orientations are called 
almost alternating orientations in the sequel, see Figure [TT1 (b). 

Proving a lower bound for the number of almost alternating Eulerian orientations yields 
a lower bound for the number of 2-orientations of G^ e . 

For the sake of simplicity we will work with alternating orientations of GjT_ 2 t-2 instead 
of almost alternating ones of In these Eulerian orientations the wrap-around edges are 
directed alternatingly up and down respectively left and right, see Figure [TT] (c). It is easy to 
see that this is a lower bound for the number of almost alternating orientations of GT„. Since 
we are interested in an asymptotic lower bound there is no difference in counting alternating 
orientations of GT_ 2 and GT« from our point of view and we will continue working with 
alternating orientations of GT« to keep then notation simple. 

Consider a vertex column Vj of G\ t and the edge columns E^_ 1 and Ej . Let X\, X2 be 
orientations of Ej_ 1 respectively Ej. Let 6(X\, X2) = 1 if and only if the edges induced by 

can be oriented such that all the vertices of have out-degree 2. Let 8u{X\, X2) = 1 
respectively 5d{X\, X2) = 1 if and only if 5(X\, X2) = 1 and the wrap-around edge induced 
by Vj is directed upwards respectively downwards. Note that 

S U {X U X 1 ) = 1 = 6 D (X 1 , X x ) 

and 

5u(X u X 2 ) = l^ 5 D (X 2 , X{) = 1. 

We define two transfer matrices Tjj(2k) and Tjj(2k). These are square 0-1-matrices with 
the rows and columns indexed by the ) orientations of an edge column of size 2k, that have k 
edges directed to the right. The transfer matrices are defined by (Tu(2k))x 1 ,x 2 = 8\j{X\, X2) 
and (T D (2k)) Xl ,x 2 = S D (Xi, X 2 ). Hence T v (2k) = T D (2k) T and T 2k = T v (2k) ■ T D (2k) is a 
real symmetric non-negative matrix with positive diagonal entries. From the combinatorial 
interpretation it can be seen that T 2 k is primitive, that is there is an integer I > 1 such that 
all entries of T| fe are positive and thus the Perron-Frobenius Theorem can be applied. Hence, 
T2k has a unique eigenvalue h.2k with largest absolute value, its eigenspace is 1-dimensional 
and the corresponding eigenvector is positive. 
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Let Xa be one of the two edge column orientations that have alternating edge directions 
and &a the vector of dimension ( 2 k ) that has all entries but the one that stands for Xa, 
which is 1. 

The number CA(2k, 2£) of alternating orientations of G 2k2e is (T 2k ) XA Xa = (eA, T 2k eA)- 
Since the eigenvector belonging to is positive it is not orthogonal to any column of T 2 k 
and we obtain ^ 

lim c A (2k,2£) 1 / e = lim ( fc) ) = A 2k , 

where the last equality is justified by an argument known as the power method. 

It follows from [22] that the limit lim^oo A 2 ( k exists, but for the sake of completeness we 
provide an argument from [TJ. We use that A 2k > (v,T 2k v) / (v,v) for any vector v and that 
(eA, T 2k eA) = (eA,T 2p eA) since both expressions count the number of alternating orientations 
°f G T 2K2p - 

4L k Y = ^ P 2 k ) 1/k > {(eA,T% k e A )) 1/k = ((e A ,T 2 k p e A )" 1/k 
Taking limits with respect to k on both sides yields 



liminf A^V > liminf ((e A , T k e A )) V * = A 2p 

k— >oo I k—*oc \ J 



1/k 1/v \lk 

which implies lim inf A 2k > lim sup A 2 . It follows that lim A 2k exists. Similar arguments 



as above yield the following. 

(e A T^T^ k e A ) _ (e A ,T^e A ) _ (e A ,T k p+Aq e A ) 



Al> 



2k ~ (Tl k e A ,n k e A ) (e A ,T%e A ) (e A ,T k q e A ) 

Taking limits with respect to k on both sides yields 

lim K y k >(^A 1,p 



) x > m = lim Al{ 4k 

k^oo 

of G 2k 2i . Using a Mathematica program we have computed Aio and As with the result that 



We are interested in lim lim c A (2k, 2£) 1/ike = lim A 2 ' k since Ak£ is the number of vertices 

k^oo t— »oo fc^oo 



^V" >( 23358714 V /4 > 1.537 
A 8 J ~ I 418.2717 J 



□ 

Remark. We return to the correspondence between 2-orientations of and Eulerian 
orientations of G\^ that was mentioned at the beginning of the last proof. By the pigeon 
hole principle, there must be a sequence of orientations X k £ of the wrap-around edges that 
extends asymptotically to (8 • \/3/9) fc ^ Eulerian orientations of G\ t This implies that for 
k,£ big enough there is an a^e on Gk/ such that there are (8 • V3/9) M otk ^-orientations of 
Gk,e- This ctk,e satisfies ctke(v) = 2 for every inner vertex v and a^iiw) G {0, 1,2} for every 
boundary vertex w. We call ce-orientations of this type inner 2-orientations of the grid. inner 2- 

We think that has asymptotically (8 • y/3/9) ki 2-orientations. But we were not able orientations 
to show this, just like for the case of the triangular grid, see the last remark of Section [3.11 of the grid 
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Theorem 9 Let Q n denote the set of all plane quadrangulations with n vertices and Z(Q) 
the set of 2- orientations of Q £ Q n . Then, for n big enough 

1.53" < max \Z(Q)\ < 1.91". 

QeQ n 

Proof. The lower bound is that from Proposition^ An upper bound of 2" follows immediately 
from Lemma [TJ Note that we may assume that Q does not have vertices of degree 2, because 
their incident edges would be rigid. Bonsma [61 [5] shows that triangle-free graphs of minimal 
degree 3 have a spanning tree T with more than n/3 leafs. As Q is bipartite, T has as a set 
/ of at least n/6 leafs which is an independent set of Q. As in Proposition [IJ this yields that 
there are at most 2 n • (3/4) n / 6 < 1.91™ 2-orientations of Q. □ 



5 Counting Bipolar Orientations 

We first give an overview of the definitions and facts about bipolar orientations that we need 
in this section. A good starting point for further reading about bipolar orientations is |llj . 

Let G be a connected graph and e = st a distinguished edge of G. An orientation X of 
the edges of G is an e-bipolar orientation of G if it is acyclic, s is the only vertex without e-bipolar 
incoming edges and t is the only vertex without outgoing edges. We call s and t the source orientation 
respectively sink of X. There are many equivalent definitions of bipolar orientations, c.f. |llj . 

Figure 15.21 shows examples of plane bipolar orientations. Note that in this case it is 
enough to have vertices s, t on the outer face, they need not be adjacent. The following 
characterization of plane bipolar orientations will be useful to keep some proofs in the sequel 
simple. 

Proposition 6 An orientation X of a planar map M with two special vertices s and t on 
the outer face is a bipolar orientation if and only if it has the following two properties. 

(1) Every vertex other than the source s and the sink t, has incoming as well as 

outgoing edges. 

(2) There is no directed facial cycle. 

Furthermore, the following stronger versions of the above properties hold for every bipolar 
orientation. 

(l'J At every vertex other than the source and the sink, the incoming and outgoing edges 
form two non-empty bundles of consecutive edges. 

(2') The boundary of every face has exactly one sink and one source, i.e. consists of two 
directed paths. 

We omit the proof that properties (1) and (2) imply that X is a bipolar orientation. The 
proof that every bipolar orientation has properties (1') and (2') (and thus properties (1) and 
(2) as well) can be found in in [40] or |35j . 

The fact that properties (1') and (2 ; ) characterize bipolar orientations of planar maps 
yields a bijection between bipolar orientations of a map M and 2-orientations of the angular 
map, i.e., of the map M on the vertex set V U T , where where T is the set of faces of M, and 
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Figure 12: Two bipolar orientations of the same graph with different out-degree sequences and 
the corresponding a-orientation of the angle graph (vertex for the unbounded face omitted). 

edges {v, /} for all incident pairs with v € V and F E T . This bijection was first described 
by Rosenstiehl [30] . 

Since bipolar orientations and 2-orientations of quadrangulations are in bijection, we ex- 
plain now how the results from Theorems l9l and [TT1 are related. A triangulation with n vertices 
has an angle graph with roughly 3n vertices. Hence the upper bound of 1.91 3n , that Theorem! 
yields for the number of bipolar orientation, is worse then the upper bound from Theorem II 11 
Conversely, every quadrangulation Q with n vertices is the angle graph of the map obtained 
by connecting two vertices of one of its partition classes by an edge, if they lie on a common 
4-face of Q. This might yield a multi- graph if Q has degree 2 vertices. But we may neglect 
this, since parallel edges must have the same direction in every bipolar orientation. One of 
the partition classes of Q has size at most n/2, and thus the upper bound from Theorem [TT] 
yields that Q has at most 3.97 n / 2 2-orientations, which is worse than the bound from Theo- 
rem The example of the G^g, which has 1.53 n 2-orientations is the angle graph of a graph 
which has roughly n/2 =: n' vertices. Therefore this yields only an example with 1.53 2n ' 
bipolar orientations, which is far away from the bound given in Theorem II 1[ Conversely, the 
triangular grid Tj^, which has at least 2.91 n bipolar orientations has an angle graph with 
roughly 3n = n' vertices. This yields a quadrangulation with 2.91™ / 3 2-orientations, which is 
worse than the bound for the number of 2-orientations that we obtained in Proposition 

5.1 Counting Bipolar Orientations on the Grid 

We now turn to analyzing the number of bipolar orientations of Gk/, with source (1, 1) and 
sink (k, i) if k is odd and sink (k, 1) if k is even. For the proof of the Theorem, we need sparse 
sequences. A sparse sequence is a 0-1-sequences without consecutive Is and it is well known sparse 
that there are F n+ 2 such sequences of length n, where F n+ 2 denotes the (n + 2)th Fibonacci sequence 
number. 

Theorem 10 Let and B{Gk/) denote the set of bipolar orientations of G^/- For k,£ big 
enough the number of bipolar orientations of the grid Gk/ is bounded by 

2.18 M < \B(G k>t )\ <2M9 ke . 

Proof. We first prove the lower bound with an argument using directed cycles in a canonical 
orientation, as in Proposition [3l To make this tool applicable we consider 2-orientations of 
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the angle graph Gk,i of Gk,i- Figure [131 shows the angle graph 64,5. The graph Gk,e has 
2k£ — 3(k + i) + 4 squares. All dotted edges are rigid, just like the four edges which are 
adjacent to a degree 2 vertex. Therefore, we may neglect all these edges in the rest of the 
proof. 

The independent set I of directed cycles in the canonical orientation is marked by dots 
and includes approximately half of all squares. The set I' consists of all squares that are not 
in /. Members of I' can be flipped if either the two cycles of I above it, or the two cycles of 
I below it are flipped, that is in 2 out of 16 cases (1 out of 4 for boundary squares). Roughly 
half of all squares are in I'. 

(b) 






Figure 13: The grid G^ is shown with its angle graph in blue; a canonical 2-orientation on 
G4.5 where red edges all connect to an additional vertex and the dots mark an independent 
set of directed cycles; the central part of G§$ and the traversal used in the proof of the upper 
bound. 

Thus, there are at least 2' / ' + ' J '' //8 bipolar orientations of Gk/, which leads to an asymptotic 
lower bound of 2 9M / 8 ra 2.18 w . 

For the proof of the upper bound we use a bijection discovered by Lieb [22]. The bijection 
relates face 3- colorings where no two squares sharing an edge have the same color and inner 

2- orientations of the square grid as shown in Figure [14] (a). Figure [14] (b) shows the face 

3- coloring corresponding to the canonical 2-orientation, that we used for the proof of the 
lower bound. 



(a) 



RiG RfB 
B|R B|G 

GfR G|B 




face 

3- colorings 



Figure 14: Lieb's bijection between inner 2-orientations and face 3-colorings on the grid. 



Here we use this relation on Gj-^. We prove an upper bound for the number of face 
3-colorings of Gf. £■ The right part of Figure [13] shows the central part of G^/ bounded by 
a thick polygon. We will encode the 3-coloring on the faces of the central part of G^i as a 
sparse sequence a, where a% represents the ith square on the path P indicated by the arrows 
in the figure. 

The set T> of faces, which are not in the central part, has less than 3' 15 ' 3-colorings. In the 
encoding described next, the code for the zth face of the path P depends only on faces in D 
and faces of P with index smaller than i. Figure [T5l shows how the color of the highlighted 
face is encoded by a or a 1. The arrows indicate the direction in which we traverse the 
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central part of the graph. There are three cases, one for a face where the path makes no turn 
and two for the two different types of turn faces. The variables X, Y, Z represent an arbitrary 
permutation of R,G, B. 

As for the decoding, it is clear from the figure that the faces marked with an X or Y plus 
the 0-1 encoding uniquely determine the color of the face in question. Thus the encoding 
is injective. It remains to show that there cannot be consecutive Is in this sequence. This 
follows from the observation, that writing a 1 means, that the two faces that will be used 
for the encoding of the next face on the path have different colors. Thus, this face will be 
encoded by a 0. 
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Figure 15: Encoding a 3-coloring by a sparse 0-1— sequence. On the left the encoding for a 
square where the path makes no turn, in the center and right for the two different kind of 
turn faces. 



We bound the number of such encodings from above. The set T> can be covered by at 
most four horizontal plus four vertical rows of faces, thus \T>\ < 4(k + £). The length of the 
path is bounded by the number of bounded faces of Gk,e, which is less than 2k£. Therefore, 
there are at most 

3< fc+£ ) • F 2ke+2 

such encodings. Using the asymptotics for the Fibonacci numbers this implies, that there are 
less than 2.619 fc ^ such encodings for k,£ big enough. □ 

Lieb's analysis of the number of Eulerian orientations of Tm is of interest in this case as 
well. It allows to improve the upper bound for grids with side lengths ratio one to two. 

Proposition 7 For k big enough the number of bipolar orientations of the grid Gk,2k * s 
bounded by 

2.18 2fc2 <|£(G fc , 2fc )|< 2.38 2fc2 . 

Proof. By G' k t we denote the graph obtained from Gk,e by deleting Voo and all incident 
edges, which are shown as dotted edges in Figure [T3j Figure [16] shows how to cut Gj 4 in 
two steps such that the grid looks like G' 3 5 (if we do not identify vertices). The last drawing 
shows, that every ci-orientation of G' 3 5 yields a Eulerian orientation of Gj 4 when we do the 
appropriate identifications. In general, this approach yields an injection from the bipolar 
orientations of Crfc+i^fe+l t° the Eulerian orientations of G\ k 2fe . As Lieb [22] has shown that 
G\ k 2k has asymptotically (8 • \/3/9) 4fe2 Eulerian orientations, this yields an upper bound of 
(64/27) 2fc2 for the number of Eulerian orientations of Gfc+i,2fc+i- Every bipolar orientation of 
Gk,2k can be complemented to a bipolar orientation of G^i 2fc+i, thus Gk2k nas a t most as 
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Figure 16: How to obtain the tilted grid £5,3 from Gj 4 with two cuts. The numbers in the 
first three drawings are the vertex labels, in the last one they indicate a. 

many bipolar orientation as Gk+i^k+i- The lower bound follows from the more general claim 
of Theorem [10l □ 

Remark. The same problems as described in the closing remarks of Sections 13.11 and 2] arise 
here when trying to show that G^^k actually has (64/27) 2fc bipolar orientations by using 
Lieb's result for the torus. 



5.2 Counting Bipolar Orientations On Planar Maps 

Note that adding edges to the faces of size at least 4 of a planar map M can only increase 
the number of bipolar orientations by Proposition [6j Thus we can restrict our considerations 
to plane inner triangulations in this section. 

Theorem 11 Let Ai n denote the set of all planar maps with n vertices and B(M) the set of 
all bipolar orientations of M £ M n - Then, for n big enough 

2.91™ < max \B(M)\ < 3.97 n . 
MeM„ 

For the proof we need a couple of facts about Fibonacci numbers, which are summarized in 
the following lemmas. The Fibonacci numbers are the integer series defined by the recursion 

F 1 = 1, F 2 = 1, F n = F n _! + F n _ 2 for n > 3. 

Define F Q = and let (j) = ±^ be the Golden Ratio. 

Lemma 3 The Fibonacci numbers have the following properties. 

V TP _ <l> n 

• llULfj—xx, r n — 

• ELo ^ p n-i = \ (n{F n+1 + F n _i) - F n ) 

The first two are standard results from the vast theory of Fibonacci numbers. The last formula 
is attributed to Shiwalkar and Deshpande in [34, A001629]. The next lemma summarizes facts 
about sparse sequences. 
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Lemma 4 The number of sparse sequences of length n is F n+ 2- Let r n (i) be the number of 
sparse sequences of length n whose ith entry is 1. Then, 

• r n {i) = Fi ■ F n+ i_.j 

n 
i=l 

. l im TJUrn(i) = 1 ^ Q 2?64 

The first identity follows from a construction of sparse sequences of length n from sparse 
sequences of length n — 1 plus the string "0" and sparse sequences of length n — 2 plus "01" . 
The second and third identity then follow using the facts from Lemma El 

Before proving the Theorem II II we give two results for the number of bipolar orientations 
of special classes of planar maps. 

Proposition 8 A stacked triangulation with n vertices has 2"~ 3 bipolar orientations. 

Proof. The K± has two bipolar orientations for fixed source and sink. We proceed by induction 
and assume, that a stacked triangulation with n vertices has 2 n_3 bipolar orientations. Now 
let T be a stacked triangulation with n + 1 vertices and v a vertex of degree 3 in T. Then, 
T — v has 2™~ 3 bipolar orientations by induction. Now stacking v into T again, there are 
exactly two ways to complete a given bipolar orientation onT — v without violating Properties 
(1) or (2) from Proposition [6l Thus, there are 2^ n+1 ^~ 3 bipolar orientations of T. □ 



Proposition 9 Let O n be the set of all outerplanar maps with n vertices. Then, 

max \B(M)\ = F„_i m 1.618 n_1 . 

Proof. We show first that there are indeed outerplanar maps with F n -\ bipolar orientations. 
Let T := T^^t be the triangular grid with two rows. We consider bipolar orientations of T 
with source (1, 1) and sink (2,£). In every such bipolar orientation the boundary edges form 
two directed paths from (1, 1) to (2,£). We start by defining the standard bipolar orientation 
Bq of G* , which is shown in Figure 1171 




Figure 17: The standard bipolar orientation on T% £. 

In Bq the vertical inner edges are directed downwards and the diagonal ones upwards. 
Now we encode any other orientation of the inner edges by a sequence (ai)j=i... ra ' of length 
n' = n — 3, where <Zj = 1 if the corresponding edge has the opposite direction as in Bq and 
a% = otherwise. The entries come in the natural left to right order in (ai)i=i... n ' ■ We show 
that all sparse sequences of length n' produce bipolar orientations. In a sparse sequence there 
are no consecutive Is, thus out of the two inner edges incident to a vertex at most one is 
reversed with respect to Bq. This guarantees that there is no directed facial 3-cycle. As all 
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vertices have an incoming and an outgoing outer edge the resulting orientation is bipolar, 
according to Proposition [6J 

It remains to show that F n _i is an upper bound for the number of bipolar orientations of 
any outerplanar map M with n vertices. We may assume that M is a plane inner triangulation. 
The proof uses induction on the number of vertices and the claim is trivial for n = 3. Now, 
let M have n + 1 vertices and let s be the source vertex. If M has a vertex x ^ s, t of degree 
2 with neighbors v, w then the direction of the edge {v,w} determines the directions of the 
edges {x,v} and {x,w}. Therefore, M has at most as many bipolar orientations as M — x, 
that is at most -F n _i many. If all vertices but s and t have degree at least three, then s and 
t have degree 2 and the vertices of every inner edge of M are separated by s and t on the 
outer cycle. This is because the interior of the boundary cycle on n + 1 vertices is partitioned 
into n — 1 triangles, and thus two of these triangles must share two edges with the boundary, 
which yields two degree 2 vertices. 

So s is incident to only two vertices v and w, and we may assume that v has degree 3 in 
M , that is the inner edge e = {v, w} is the only inner edge incident to v. Now, let X be some 
bipolar orientation of M in which e is directed from v to w. Then, the orientation of M — s 
induced by X is a bipolar orientation with source v. For a bipolar orientation Y in which e 
is oriented from w to v, the orientation of M — s induced by X is a bipolar orientation with 
source w and v is a vertex of degree 2 in M — s. This mapping is injective, and thus M has at 
most as many bipolar orientations as M — s and M—{s, v} together, that is F n _i+F n _2 = F n . 

□ 

Remark. From the above proof it also follows, that T^i is the only outerplanar map on 21 
vertices, which has bipolar orientations. 

The example that gives the lower bound for the number of bipolar orientations of planar 
maps is the triangular grid % with source (1, 1) and sink (k, k). 

Proposition 10 Let T^k be the triangular grid and k big enough. Then, 

\B(T htk )\ > 2.91". 

Proof. We first claim that T^j. has at least 2.618 fc2 bipolar orientations. To see this we glue 
together k — 1 copies of T^^,- Every orientation of T^k obtained in this way corresponds to 
a concatenation of k — 1 sparse sequences of length 2k — 3, which we call an almost sparse 
sequence. We denote the set of all such sequences of length 2k 2 — 5k + 3 by S, the cardinality 
of S is -F^fc-i wri ich is bounded below by i ? 2fc 2 -5fc+3 — 2.618 fc for k big enough. That each 
s £ S corresponds to a bipolar orientation of T^k can De checked using Proposition [6l 

The horizontal edge e$ ,- := — ► + 1) lies on the boundary of two triangles for 
2 < i < k — 1. The other four edges of these triangles are 

{(»,j),(i-l,j + l)}, + + + + + 

The crucial observation for improving the above bound is, that we can reorient ejj if and only 
if the entries belonging to these four edges show one of the two patterns 10 ... 01 or 01 . . . 10. 

We now choose k — 1 sparse sequences of length 2k — 3 independently uniformly at random 
and concatenate them to obtain a random almost sparse sequence s G S. It follows from the 
first identity from Lemma U] with n = 2k — 3 and i = 2j — 1 that for + there are 

^2j-l-^2fc-2j-l-^2fc^o. sec l uences that have {(i,j), (i — 1, j + 1)} marked 1 out of the total F^k-i 
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sequences. This is used to calculate the probability that the entry for — + 1)} is 

1 as 



Ms) = i]>^ 



Hm F 2 j _ i F 2 fc 2 j - 1 ^2k-i = Um 1 . 2i-l . . -2fc+l = 1 (5) 

Taking the limit is only justifiable if (2j— 1) — > oo and (2/c— 2j— 1) — > oo for — > oo. Therefore 
we introduce 5 > and denote the set of horizontal edges with S(k — 1) < j < (1 — S)(k — 1) 
and 2 < i < k — 1 by i£<j. Taking the limit in equation (5) is justified for all e$ G The 
size of this set is \E 5 \ = (1 - 25) (fc - l) 2 . 

This works analogously for the edge {(i, j + 1), (i + and the events are independent. 
The pattern 01 ... 10 has the same probability in the limit and the patterns mutually exclude 
each other. Thus, for every e > the probability that the edge + 1)} € Eg can be 

flipped is 

2 

for k big enough. 

The rest of the proof is to analyze how many of the flip-patterns we expect for a sparse 
sequences s. Let Q(s) = Yli,j ^-i,j( s ) be a random variable counting the number of flippable 
edges in s. We use Jensen's inequality to estimate the number of orientations, which is at 
least 

\S\ ■ E S&S [2W] > F 2k2 _ 5k+3 ■ 2 E ^[Q( S )1 > 2.618 fc2 • 2 (1 ^ ){ ^ 1)2 (^H > 2.91 fc2 

for k big enough. □ 

Remark. The third identity from Lemma U] gives the expected number of Is that a random 
entry of a random sparse sequence of length n has. In equation ([5]) we calculate the expected 
number of Is in a random sequence of length n at a fixed entry which is far enough away from 
the boundary. The fact that these two values agree shows that the expected number of Is at 
a fixed entry of a sparse sequence does not depend strongly on the choice of the entry. 

The following relation is useful to upper bound the number of bipolar orientations for 
plane inner triangulations. It has been presented with a different proof in [25J. Let J-f, be 
the set of bounded faces of M and B the set of bipolar orientations of M. Fix a bipolar 
orientation B. The boundary of every triangle A € consists of a path of length two and a 
direct edge from the source to the sink of A. We say that A is a + triangle of B if looking 
along the direct source-sink edge the triangle is on the left. Otherwise, if the triangle is on 
the right of the edge we speak of a — triangle, see Figure [T8l We use this notation to define 
a mapping G B ■ F b -> {-, +}. 

Theorem 12 Let M be a plane inner triangulation and B a bipolar orientation of Ad. Given 
Gb, i.e., the signs of bounded faces, it is possible to recover B. In other words the function 
B — ► Gb is injective from B(M) —> { — , -j-jl^l _ 

Proof. Given we construct B. We start by orienting all edges on the boundary of the outer 
face such that s and t are the unique source and sink of this face. We extend this partial 
orientation Y with two rules. The vertex rule is applied to a vertex v, which already has 
incoming and outgoing edges. It takes a bundle of consecutive edges of v which is bounded 
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Figure 18: A bipolar orientation, the corresponding +/— encoding and an illustration of the 
decoding algorithm. 

by two outgoing edges. It orients all the edges of the bundle such that they are outgoing 
at v. Note that these edge orientations are forced by Property (1') of bipolar orientations in 
Proposition El The face rule is applied to a facial triangle A, which has two oriented edges. 
The sign Gb(A) is used to deduce the orientation of the third edge. 

Note that these two rules preserve the property that every vertex v, which is incident to 
an oriented edge in Y, can be reached from s along an oriented path. In particular v has an 
incoming edge. 

Let Ay be the union of all faces, which have all boundary edges oriented. Initially, 
Ay consists of the outer face. Since B is acyclic the boundary of Ay is acyclic as well. 
Consequently, as long as there are faces which do not belong to Ay, there is a vertex v on the 
boundary of Ay which has two outgoing edges that belong to the boundary of Ay. Either v 
is a candidate to extend the orientation using the vertex rule or there is a face incident to v 
which becomes an element of Ay by applying the face rule to it. 

We have thus shown that the rules can be applied until Ay is the whole plane, i.e., all 
edges are oriented. They have to be oriented as in B, by construction. □ 

The next theorem gives a necessary and sufficient condition for a vector in {— j+j'^' to 
induce a bipolar orientation. For the sake of simplicity we state it only for triangulations, 
but the generalization to inner triangulations is straight forward. In order to obtain a more 
elegant formulation, we adopt the convention, that the unbounded face is signed + if the 
bounded face adjacent to sink and source is signed — . Otherwise the unbounded face is 
signed — . Thus we work now with signings of the set T of all faces. We say that a + triangle 
is the right knee of the vertex at which it has in incoming and an outgoing edge. Similarly a 
— triangle is the left knee of exactly one of its vertices. 

Theorem 13 Let T be a triangulation, x & {— , and T~ and J- + the sets of faces that 

have negative respectively positive sign in x. Let M(T) + and M(T)~ denote the subgraphs of 
the reduced angle graph M(T) — {s, t} induced by V U T~ respectively V U J- + . 

Then, x induces a bipolar orientation on T if and only if both, M(T) + and M(T)~ , have 
a unique perfect matching. 

Theorem [13] implies that every vertex of T other than s,t must be adjacent to at least 
one + triangle and one — triangle. 

Proposition 11 Let T be a plane inner triangulation. Then, T has at most 3.97 n bipolar 
orientations. 
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Proof. By Euler's formula, there are 2^~ 1 < 2 2n_2 ~^°° = 4 n ~ 1 • 2~^°° many binary vectors of 
length / — 1. By the bijection from Theorem 1121 4 n ~ 1 • 2 - -' 00 is also an upper bound for the 
number of bipolar orientations of T. 

To squeeze the bound below 4 n we use the above observation that every vertex of T must 
be adjacent to at least one plus + triangle and at least one — triangle. Thus, out of the 2 dv 
possible +/— vectors at a vertex v at least two are not feasible. Similarly, at an outer vertex 
v 7^ s,t, there is exactly one angle forming a knee at v. The sign of this angle depends on 
which of the two oriented paths of the outer boundary v lies on, but it is fixed either way. 
Thus, out of the 2 dv ~ l possible sign patterns at v at least one is not feasible. We summarize, 
that at most a fraction of (1 — 2 1 ~ dv ) of all sign vectors is potentially feasible at every vertex 
but s and t. We denote the set V \ {s, t} by V' and its cardinality by n', i.e. n' = n — 2. 

We apply Jensen's inequality, which says that for a convex function <p the inequality 
ip Xi/n) < ^{xi}) jn holds and the inequality is reversed if cp is a concave function. As 
logx is concave we obtain 

log ( ( n a - 2i ^)) 1/n ) = \, e (i - ^) ± ( i e (i - ^ 

By the monotonicity of the logarithm this implies 



n(i-2^)<(i£(i- 2 i- 



veV 

The function 2 X is convex and we apply Jensen's inequality again which yields 

1 E (l - 2 1 ~ £ M = 1 - 2 ( — t 2 ~ dv ) < 1 " 2 1+l l n '^v 
n v<=v \ veV J 

Since we deal with simple plane graphs YlveV ^ — 2(3n — 6) = 6n' and we conclude 

n a - 2^) ^(i-^ 6 )-' =(§)"'. 

Since we deal with plane graphs, T can be partitioned into at most four independent sets 1^, 
k = 1, ... ,4 by the Four Color Theorem. Thus, 

4 



I)" ^ iK 121 d )=n IK 1 - 21 

' ' veV k=ivei k 



and for at least one of the independent set it must hold that 

n'/A 

We are ready to conclude that there are at most 

/oi\(n-2)/4 /S?\ 1/2 
4 n-l . 2 -/oc < 3 Q7 n . 2 -/oc l±_\ < 3 9? n 

bipolar orientations of T. □ 
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6 Complexity of Counting a- Orient at ions 



Given a planar map M and some a : V —* N, what is the complexity of computing the number 
of a-orientations of Ml In some instances this number can be computed efficiently, e.g. for 
perfect matchings and spanning trees of general planar maps. In Section 16.11 we show that 
counting is #P-complete for other a-orientations, and in Section T6. 21 we discuss how to adapt 
an existing fully polynomial randomized approximation scheme. 

6.1 #F-Completeness 

Recently, Creed [8] has shown the following theorem. As already mentioned Eulerian orien- 
tations are a-orientations and hence Theorem Q3] says that counting a-orientations is H=P- 
complete. Since we use the proof technique in the sequel, we sketch the proof here. It uses a 
reduction from counting Eulerian orientations, which has been proven to be #P-complete by 
Mihail and Winkler in |26j . 

Theorem 14 It is #P '-complete to count Eulerian orientations of planar graphs. 

Proof. We aim to show that number Eulerian orientations of a graph G can be computed 
in polynomial time with the aid of polynomially many calls to an oracle for the number of 
Eulerian orientations of a planar graph. 

In order to count the Eulerian orientations of a graph G with n vertices a drawing of this 
graph in the plane with I crossings is produced. We may assume that no three edges cross in 
the same point and that I is of order 0(n 2 ). 

From this drawing a family of graphs Gi for i = 0, . . . , £ is produced. In Gi every crossing 
of two edges {u, v} and {x, y] is replaced by the crossover box Hi on 4k + 1 vertices, see 
Figure QjJl For example in Go every crossing is replaced by just one vertex. We call the edges 
{v,w^}, {u, w^}, {x, w%}, {y, w^} the connection edges of Hk- connection 

Every Eulerian orientation of G induces a Eulerian orientation of Go , but there are Eule- edges of Hk 
rian orientations of Go, which do not come from a Eulerian orientation of G. Given a Eulerian 
orientation of Go we call the orientation of the edges incident to a vertex wq, which replaces 
a crossing, valid if exactly one of the edges {v,wo}, {u,wo} is directed away from wo, and 
invalid otherwise. 

Hq Hi Hk 




X u x u 



Figure 19: The crossover boxes defined in [8]. 

We call a configuration of the connection edges of H^ valid if exactly one of the edges 
{v,w1}, {u,w^} is directed towards {u, v} and exactly one of the edges {x,wf,}, {y,wf.} 
is directed towards {x,y}. We call a configuration of the connection edges of Hk invalid if 
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exactly two of the connection edges are directed towards {u,v,x,y}, but it is not a valid 
configuration. By xt respectively y& we denote the number of ways that a valid respectively 
invalid configuration of the connections edges of Hk can be extended to an orientation of the 
edges of Hk, such that all vertices of Hk have out-degree 2. 

It is clear that xq = 1 = j/o an d Creed observes the following recursion 

x k = 4x fc _i + 2y k -i 
y k = 4a:fc_i + 3y fc _i. 

This recursion formula can be verified by a simple enumeration. In [8] a lemma from [39] 
is used to argue that the sequence x k / y k is non-repeating. In Lemma [5] we provide an easy 
argument to show from first principles that Xk/y k is strictly monotonically decreasing. 

Let Ni denote the number of Eulerian orientations of Go which have i valid cross over 
boxes, that is N( is the number of Eulerian orientations of G. The number of Eulerian 
orientations of Gk is 

I 

EO{G k ) = Y J N ^lVk l - 

i=0 

Hence the number EO(Gk)/yl * s * ne vame °f the polynomial p(z) = Y^i=oNiZ l at the point 
Xk/yk- Since computing yjr. is easy the polynomial p of degree I can be evaluated at I + 1 
different points with £+1 calls to an oracle for counting Eulerian orientations of planar graphs. 
Hence the coefficients of p can be determined using polynomial interpolation. In particular 
this yields a way to compute iVg. □ 



Lemma 5 Let xq = 1 = y$ and Xk+i = 4x/c + 2yk while y k +i = 4xfc + 3yk- Then the sequence 
x k /yk is strictly monotonically decreasing. 

Proof. It follows directly from the recursion that Xk+i = yk+i — Vk and yk+i = — tyk-i- 
We want to show that 

Xk+i 4x fc + 2y k Xk_ ^ / Xk \ 2 _ Xk_ _ 2 Q 



Uk+i 4x fc + 3y k y k \y k J y k 

Since Xk/yk > this inequality is satisfied if and only if Xk/y k > (1 + \/33)/8 =: c. Note that 
1/(1 — c) = 3 + 4c can be easily derived since c solves 4i 2 — t — 2 = 0. 

It remains to show that Xk/yk = 1 — (Vk-i/Vk) > c i which we do by induction. Since 
xo/yo = 1 and x±/yi > 6/7 the induction base is trivial and the induction hypothesis guar- 
antees that 1 — c > yk-2/yk-i f° r k > 2. We have 

Xk Vk-i ^ ^ Vk ^ 1 
1 > c <J=^ > 



yk yk yk~\ l - c 

Using the induction hypothesis for the last inequality we obtain that 

y k 7y k -i - 4y fc _ 2 yfc-2 ^ „ . . 

7 — 4 > 3 + 4c 



2/fc-i Vk-i Vk-i 1 - c 

□ 

The reduction from counting perfect matchings in bipartite graphs in [26j to counting 
Eulerian orientations creates vertex degrees, which grow linearly with the number of vertices 
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of the graph. This implies also that the image of a(v) = d(v)/2 grows with n. The following 
theorem shows that this dependence of a on n is not necessary to obtain a #P-completeness 
result. 

Theorem 15 For the following instances of graph classes and out-degree functions a the 
counting of a- orientations is ^P-complete. 

1. Planar maps with d{v) = 4 and a{v) E {1, 2, 3} for all v E V. 

2. Planar maps with d{v) E {3, 4, 5} and a{v) = 2 for all v E V . 

The proof uses the planarization method from the proof of Theorem 1141 in conjunction with 
the following theorem from [9]. 

Theorem 16 Counting perfect matchings of k-regular bipartite graphs is j^P-complete for 
every k > 3. 

Proof of Theorem 1 1 51 Perfect matchings of a bipartite graph G with vertex set V = AUB are 
in bijection with a-orientations of G with a(v) = 1 for v E A and a(v) = d(v) — 1 for v E B. 
This bijection is established by identifying matchings edges with edges directed from A to 
B. Hence in /c-regular bipartite graphs counting perfect matchings is equivalent to counting 
what we call l-(k — l)-orientations in the sequel. 

We observe that the planarization method from the proof of Theorem [TH can be used in 
a more general setting. Let Qd be the set of all graphs with vertex degrees inflcN and Vd 
the set of all planar graphs with degrees in D. Let IcN and associate with every G E Qd & n 
out-degree function «c whose image is contained in /. Then, the proof of Theorem [T4l shows 
that counting the a^-orientations of graphs in Qd can be reduced to counting a^-orientations 
of the graphs in Pdu{4} where the image of a' G is contained in / U {2} for all G' E Vdu{4}- 

When we apply this to 4-regular bipartite graphs with 1-3-orientations, that is perfect 
matchings, it yields the first claim of the theorem since counting perfect matchings of bipartite 
4-regular graphs is #P-complete by Theorem [161 We give two different proofs for the second 
claim. 




Figure 20: Gadgets to translate (a) a{y ) = 1 to a = 2 for d(v) = 3, (b) a(v) = 1 to a = 2 for 
d(v ) = 4, and (c) a(v) = 3 to a = 2 for d(v) = 4. 

Let G be a graph with a degree 3 vertex v and etc an associated out-degree function with 
<xg{v) = 1- We substitute v by the gadget G\ from Figure [20l (a) to obtain a graph G' . The 
gadget has five vertices that induce nine edges and 3 edges connect it with the neighbors 
of v in G. Let aa> = c^g on V{G) — v and be q.g'{u) = 2 for u E V{G\). Exactly one of 
the connection edges must be directed away from G\ in every agz-orientation of G' . Note 
that G\ is symmetric in the three connection vertices. It is easy to check that G\ has ten 
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orientations with out-degree 2 at every vertex once the outgoing connection edge has been 
chosen. Thus, every a^-orientation is associated with ten o^-orientations of G' and since the 
underlying ag-orientations can be reconstructed from every olqi -orientation we obtain that 
r aal (G') = 10 -r aG (G). 

Let a 3-regular bipartite graph G with vertex set V = A U B be given and G' be obtained 
from G by substituting every vertex v G A by a copy of G\. The number of 1-2-orientations of 
G is lO"'" 4 ' -r a ,{G'). Using the planarization method this yields that counting 2-orientations 
of graphs from Vis,A,5} 1S #-P-complete. 

Similarly counting perfect matchings of 4-regular bipartite graphs can be reduced to count- 
ing 2-orientations of graphs from Vs^^sy. We mention this second proof since it uses planar 
gadgets. More precisely G' is obtained from G by substituting the vertices of one partition 
class of a 4-regular bipartite graph G by the gadget G2 shown in Figure [20] (b) and the ver- 
tices from the other partition class by the gadget G3 shown in Figure [20] (c) . Using a similar 
reasoning as above one obtains that the number of 1-3-orientations of G is (26 -6)~^r a , (G'), 
where G2 and G3 give a blow up factor 26 respectively 6. □ 

Having proved Theorem [15] it is natural to ask whether it is #P-complete to count a- 
orientations for A:-regular planar graphs and constant a. This setting implies that a = k/2. 
Planar graphs have average degree less than six. Furthermore a 2-regular connected graph, 
that is a cycle, has two Eulerian orientations. Hence the above question should be asked for 
k = 4. Note, that the planarization method yields vertices with d(v) = 4 and a(v) = 2, so 
we do not need to restrict our considerations to planar graphs when trying to answer this 
question. 

Problem. Is it ^P-complete to count Eulerian orientations of 4-regular graphs? 
To the best of our knowledge even the following problem is also open. 

Problem. Is it #P-complete to count Eulerian orientations of graphs with degrees in 
{1, . . . , k}, for some arbitrary but fixed k £ N? 

We present one more #P-completeness result since it has a nice connection with the 
first problem stated above. Note, that for 4-regular bipartite graphs 2-factors and Eulerian 
orientations are in bijection. For 3-regular bipartite graphs 2-factors are in bijection with 
their complements the 1-factors, i.e. perfect matchings. Hence it is ^P-complete to count 
2-factors of 3-regular bipartite graphs. The next theorem generalizes this observation. 

Theorem 17 For every i > 3, i 7^ 4 counting 2-factors of i-regular bipartite graphs is #P- 
complete. 

Proof. The case i = 3 follows from Theorem [16] as explained above. The proof for i > 5 is a 
reduction from counting 2-factors of 3-regular bipartite graphs. The method comes from the 
proof of Theorem [16] in [9] . 

The following preliminary considerations will be needed later. We fix some edge eo of 
the complete bipartite graph on 2i vertices Ka. This graph has i 2 edges and every 2-factor 
of Kii has 2i edges. Let q be the number of 2-factors of Ki^. We want to find the ratio 
between 2-factors of Kii which contain eo and 2-factors which do not contain eo- Consider 
pairs of 2-factors and edges (F,e). Obviously, there are i 2 ■ Ci such pairs and we have that 



32 



\{(F, e) \ee F}\=2i-Ci while \{(F, e) \ e F}\ = {i 2 - 2i) ■ q. It is obvious that 

\E\ 

{(F,e)\eeF}=[j{(F,e j )\e j EF} 

3=1 

and symmetry implies, that all sets {(F, e,) | ej G i 7 } have the same cardinality. We conclude, 
that there are aj = 2cj/z 2-factors including eo- It follows similarly, that there are bi = 
(1 — 2/i)cj 2-factors not including eo- We infer that Oj/6i = 2/(i — 2) < 2/3 for i > 5. 

The gadget Pi (A;) is a concatenation of disjoint copies of K{{ — eo, with A; — 1 connection 
edges as shown in Figure [2T] for i = 5 and k = 4. The gadget is connected to the rest of the 
graph via two edges at the degree i — 1 vertices. 



Figure 21: The bridge gadget Ps(4) 

Let G be a 3-regular bipartite graph with 2n vertices and G'(k) be obtained from G by 
augmenting it with n(i — 3) disjoint bridge gadgets Pi(k) such that G'(k) is i-regular. Let P 
be some fixed bridge. Note that every 2- factor of G'{k) includes all or none of the connection 
edges of P. This is because a 2-factor is a partition into cycles and therefore intersects every 
edge cut of cardinality two in either none or both of the edges. 

Let cq be the number of 2-factors of G. Among the 2-factors of G'(k) let S denote the set 
of those, which induce a 2-factor of G, and S c the other 2-factors. We have |5| = cq ■ b\ ^ 
since the 2-factors in 5 cannot include any connection edges and thus can be partitioned into 
a 2-factor of G which is augmented by one of the b\ possible 2-factors in every bridge. 

Next we want an upper bound for \S C \. Every 2-factor in S c includes the connection edges 
of at least one bridge P. Since 2 3n is the number of subsets of E(G) each of which can be 
augmented to a 2-factor in at most ways on P we have that \S C \ < a^b\ 3 ^ 1 ^2 3n . The 
number of 2-factors of G'(k) is cq/^\ = \S\ + \S C \ for which we have the following bounds. 

,k-n(i—3) , , ik-nti— 3) . hi k-(n(i— 3) — l)r,3r? 

^ cg < c Gf(k) b~ k ^- 3) < + 

Since ai/bi < 2/3 we have that 

2 3n (j^j < \ for k > 3/(log 2 3 - l)n + l/(log 2 3 - 1). 

Note that the lower bound for k is linear in n. 

We finally conclude that cg = L c G'&i k \ f° r k large enough. Since b i kn ^ 3 -* is easy 
to compute and G k has size polynomial in n this proves the theorem. □ 
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Remark. We would like to point out, that the missing case k = 4 in Theorem 1171 cannot be 
fixed by substituting the bridge gadget by another gadget. It is crucial that 04 /bi < 1 and 
symmetry implies, that in every 4-regular graph the number of 2-factors including a fixed 
edge eo is equal to the number of 2-factors not containing eo- 

6.2 Approximation 

Let / : V — > N be defined on the vertex set of a graph G. An /-factor of a graph G is a 
subgraph of G such that every vertex has degree f(v) in F. 

Counting a-orientations can be reduced to counting /-factors in bipartite planar graphs 
and to counting perfect matchings in bipartite graphs. We next describe these transforma- 
tions. They are useful because bipartite perfect matchings have been the subject of extensive 
research (for example [24 1 |29 ] [20] ) . 

First, note that the a-orientations of M are in bijection with the the a'-orientations of the 
bipartite planar map M' obtained from M by subdividing every edge once. Here a'{v) = a(v) 
for the original vertices of M and a'(v) = 1 for all subdivision vertices. The a'-orientations 
of M 1 are in bijection with the /-factors of M' where f(v) = a'(v) for all vertices of M' . The 
bijection works by identifying factor edges with edges directed from a vertex of M to an edge 
vertex. 

The idea for the next transformation is due to Tutte [37] . The graph M' gets blown up 
into a graph M" such that M" has Y\ v& m(^( v ) ~ f( v )V- times as many perfect matchings, 
i.e., 1-factors, as there are /-factors of M'. To obtain M" from M' substitute v S V(M) by 
a Km v \m v -\ -f( v ) > such that each of the d{y ) edges incident to v in M' connects to one of the 
vertices from the partition class of cardinality d(v). 

In [20] Jerrum, Sinclair, and Vigoda give a fully polynomial randomized approximation 
scheme for counting perfect matchings of bipartite graphs. Thus, the above transformation 
yields a fully polynomial randomized approximation scheme for a-orientations as well. 

The number of perfect matchings of a bipartite graph with a Pfaffian orientation can be 
computed in polynomial time. Little [23J gave a full characterization of graphs with a Pfaffian 
orientation and in [29j a polynomial time algorithm to test whether a given graph is Pfaffian 
is introduced. As a special case of Little's characterizations it follows that all planar graphs 
are Pfaffian. Hence, in all cases where M" is planar the counting is easy. For spanning 
trees the above transformations yields planar graphs, while for Eulerian orientations it does 
not, as Theorem 1141 implies. Although we do not have a hardness result for Schnyder woods 
and bipolar orientations on planar maps, there are in both cases instances for which the 
transformation yields a non-Pfaffian graph. One such example is the augmented triangular 
grid from Section 13.11 Figure [22] shows a local structure with five vertices and four faces that 
implies that M" is not Pfaffian. The figure shows simplified versions of M' and M" . We can 
choose M' = M since the angle graph is bipartite. The Tutte transformation substitutes face 
vertices by and primal vertices of degree d by K 2j d- Instead one can simply create a copy 
of every primal vertex with the same neighborhood as the original and leave the face vertices 
unchanged to obtain a simplified version of M" . 
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Figure 22: An oriented subgraph of a triangulation with a bipolar orientation that induces a 
central 1^3-subdivision. 

7 Conclusions 

In this paper we have studied the maximum number of a-orientations for different classes of 
planar maps and different a. In most cases we have exponential upper and lower bounds 
and c|} for this number. 

The obvious problem is to improve on the constants cl and cjj for the different instances. 
We think, that in particular improving the upper bound of 8 n for the number of Schnyder 
woods on 3-connected planar maps is worth further efforts. For bipolar orientations the more 
efficient encoding from Theorem 1121 helps to improve the upper bound. We think that finding 
a more efficient encoding for Schnyder woods might be needed to substantially improve on 
the 8 n bound. 

Results by Lieb [22j and Baxter [1] yield the exact asymptotic behavior of the number 
of Eulerian orientations for the square and triangular grid on the torus. This yields upper 
bounds for the number of 2-orientations on the square grid and the number of Schnyder woods 
on the triangular grid. We have not been able to take advantage of these results for improving 
the lower bounds for the number of 2-orientations respectively Schnyder woods. 

For some instance of a-orientations there are #P-completeness results. This contrasts 
with spanning trees and planar bipartite perfect matchings for which polynomial algorithms 
are available. It remains open to determine the complexity of counting Schnyder woods and 
bipolar orientations on planar maps and of counting Eulerian orientations of graphs with 
bounded maximum degree. 
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